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Abstract 

The phase transition of Reissner-Nordstrom AdS^ interacting with a massive charged scalar field 
has been further revisited. We found exactly one stable and one unstable quasinormal mode regions 
for the scalar field. The two of them are separated by the first marginally stable solution. 
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The AdS / CFT correspondence [l[ has been used as a very important tool to probe physics 
beyond the usual application in stringtheory, namely quark gluon plasma, fluid mechanics 
as well as condensed matter models [2|-l5|. In the present discussion we analyse a model 



first considered in 



4J and later by many authors |5Ml(| related to superconductors and 
corresponding phase transitions, generally referred to as holographic superconductors. 

The AdS / CFT correspondence relates the gravitational fields in a bulk to properties of a 
field theory at its border. Moreover, the weak coupling in one of the counterparts becomes 
the strong coupling limit of the other. In {4, 11] a break of the U(l) gauge symmetry has 
been achieved by a charged scalar interacting with a black hole in the bulk, a mechanism 
leading, eventually, to a superconductor at the boundary — what has been called a holo- 
graphic superconductor. Several similar mechanisms have been proposed in the literature, 
n particular, it has been shown that several phase transitions may occur in a given model 
4j. Such a conclusion has been formulated upon computing the perturbations in the bulk 
and looking for the time independent perturbations, which presumably divide sectors of 
growing and decaying modes. 

Here we show that there are just two sectors, namely a superconducting and a normal 
phase, in view of the fact that even beyond the marginally stable solution growing modes 
still exist. Such a result comes about from the existence of further unstable modes, which 
we were able to compute explicitly by the Quasi Normal Modes technique 

We consider the bulk action of a massive charged scalar field 1/) interacting with an Abelian 
field generated by a charged AdS black hole. The Lagrange density is that of Einstein 



12|. 



gravity with a cosmological constant, a charged scalar and a gauge field, that is, 
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The above Lagrangian has also been used in (4] in the same context. The background 
spacetime is a spherically symmetric Reissner-Nordstrom AdS 4, that is, 



with 



ds 2 = -fir) dt 2 + -^-dr 2 + r 2 dtll 
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where Q and M are respectively the electric charge and mass of the black hole. Besides, L 



refers to the AdS radius, which is related to the cosmological constant A by L = a/— 3/A. 



The Hawking temperature of this black hole is given by 
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One identifies the parameters of the black hole with a critical temperature T* that charac- 
terize the phase transition. In this T* can be related to the critical temperature T c of 
a holographic superconductor. 

The solution for the gauge field reads 



A^dx 11 = $>(r)dt 



(5) 



where the electric potential is written in terms of black hole charge Q and the position of 
event horizon r + as 

• =«-« . (6) 
r r + 

For our purposes, we assume that the scalar field ip is a small perturbation that does not 
backreact. In this regime, we can use the quasinormal modes techniques in order to compute 
the "0-propagation. equation of motion results in 



D^D^ifj = m 2 ip 



(7) 



where = — iqA^. The parameters m and q are the mass and charge of the scalar 
field, respectively. Expanding Eq.flTJ) we find 



Utj) - 2iqA fl g^d v ^) - q 2 A fl A 1 *^ - m 2 vb = 



(8) 



The s-wave mode (1 = mode in the angular momentum expansion) is easily obtained. 
For our purposes it is sufficient since higher angular momentum modes are connected with 
higher decay modes, and we only need the firs^unstable mode to install unstability, as we 
are going to argue more generally. We follow 
are m, q and L. We find 



4| and set r + = Q = 1 and our parameters 
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We use the tortoise coordinate dr* 
equation 
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y and the Ansatz ip = * for the field, resulting in the 



8t 2 



+ 2iq$ 



dr 2 "* dt 



V(r)W = 



(10) 



with 



V(r) = f(r) 
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An analysis of Eg. ( ITU]) shows that the procedure used in [uj| to discretize the spacetime 
and integrate the equation of motion is not convenient because of the presence of the term 
2iq&^-. The method used in the mentioned work normally has numerical error proportional 
to A 4 , where A is the step used, but the extra term has an error proportional to A 2 , so 
we used the finite difference method, whose error is proportional to A 2 without the need of 
another coordinate change. 

Finite difference method 

We define V(r*,t) = V(-jAr*,lAt) = V jjh V(r(n)) = V(-jAn) = Vj and $(r(r*)) = 
$(— jAr*) = $j to rewrite equation (fTUj) as 
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With a Gaussian distribution with finite support as initial condition and Dirichlet conditions 
at r* = 0, we derive the evolution of \l/ by 
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The Von Neumann stability conditions usually require that < 1. If the effective 
potential Eq. (fTTj) is too large, the method is unstable even for small At. In a simple 
wave equation x + w 2 x = 0, this requirement means that the step At must be smaller than 
uj~ 1 . In our case, V(r) is proportional to jz (2 + m 2 L 2 ) for large r (i.e. r* close to 0), and 
for m 2 L 2 = 4 we can see that even with ratios as low as = 0.7 the method becomes 
unstable, so we used = 0.5. Therefore, in this work the stability conditions depend on 
the parameters m and L. 

Horowitz- Hubeny method 

Another method to find quasinormal modes in asymptotically AdS spacetimes was de- 
veloped by Horowitz and Hubeny [l_5| and used in [la, Il7 |. 



Rewriting the metric Eq.(j2]) and the gauge field Eq.([3]) in terms of ingoing Eddington- 
Finkelstein coordinate, v = t + r#, results in 



ds z = -f(r) dv 2 + 2 dv dr + r 2 dVL 2 2 , A v dv = §(r)dv 



(14) 



In the new coordinates and after performing the separation of variables ip(v,r,6,(j)) 
Yi m (9,(p) e~ luv the wave equation Eq.([9]) reads 



f^ + \f>-2i{ U + q*)] f-V(r)Z = 



where the potential is 



V(r) = [m z + 
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Following the procedure outlined by Horowitz and Hubeny 
the wave equation ([T5]) as 
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where we were be able to put the functions s(x), t(x) and u(x) in a polynomial form 
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where we set I = and Q — 1. 

Expanding the solution Z(x) to the wave equation ffTTj) around the event horizon x + = 
l/r + , we have 

oo 

Z[x) = ^a n {u){x-x + ) n . (23) 
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(18) 
(19) 
(20) 

(21) 
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The key point is to compute de roots of the equation Z(x = 0) = following from the 
boundary condition at spatial infinity. Actually, we have to truncate the de sum at an 
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intermediate n = N and check that for the greater n in which the root converge. We 
compute the zeros uj of Y^nLo a n{u)(— x + ) n using the software Mathematica and also the 



routine zroots 
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Discussion and Results 
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Figure 1: Behavior of effective potential for several values of L (color online). 



In Figure (CD) we plot the effective potential fTTTT) due to the charged scalar field propagation 
in the geometry ([2]) since the stable region, passing through the transition point next to 
L = 1.28 into the unstable region. We clearly see that the potential defines a positive 
definite barrier near the event horizon (r + = 1) and becomes negative at intermediate 
values of radial coordinate diverging at spatial infinite as expected for asymptotically AdS 



spacetimes. Since we have an additional term in Eq. fTTOj) . the result from [15| that a positive 
definite potential implies stability is not necessarily valid. However, modes with positive 
definite potentials are stable in this work. Besides, the first unstable mode is presumably 
related to the first value of L (or highest value of T) whose potential admits a bound state 
with negative energy. 

Figure ([2]) shows five different behaviors of %p(r*,t) given by Eq. (fl3|) at r* = —5 for 
m 2 L 2 = 4 and qL = 10. The results show an interesting phenomena. At a value next to 
L = 1.28 we have a phase transition. This is the same transition obtained by Gubser. Our 
results agree with those, that is, beyond this value of L the black hole become unstable. 

Fitting \1/ = Aexp(ujit) cos(u r t + 5), we get the quasinormal frequencies. We cannot tell 
the sign of u r using this fitting, but if we assume that u r depends smoothly on L, a sign 
change would appear as a discontinuity in the first derivative, as our data appears. Choosing 
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Figure 2: Quasinormal modes versus t (color online). 



the sign of u r as the same of cuj, we obtain the figure (j3J), which shows that both cu r and cjj 
change sign at the same value of L. In the same figure we adjust the scale so that the sign 
change is clearly seen. We are also interested in the dependence of these frequencies with 
the temperature, shown in figure (j4j), in this figure we also adjust the scale to better see 
the transition. We noticed that the transition occurs at a value of L smaller than expected 
in {4]. Since the numerical method used to solve second order partial differential equations 
has errors larger than methods commonly used to solve second order ordinary differential 
equations, we believe that this difference is due to numerical errors. 

There are further secondary quasinormal modes. We found the next one which is pre- 
sented in figure (jSJ). Indeed, the secondary mode shows a zero frequency at L next to the 
second marginally stable mode found in |4j. However, the first mode remains, and it is 
unstable. This means that the model remains unstable and superconductivity still holds 
beyond that second marginally stable mode. 
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Figure 3: First mode frequencies (left) and near the transition (right), (color online) 




The frequencies of the second mode displays a behaviour similar to the first mode. Again, 
choosing the sign of u r as the same of Ui, we see that both frequencies change sign at the 
same point. These frequencies are shown in Figure fl6]). 

In Figure ((Tj) we have the comparison between the results obtained using two different 
numerical approaches: the finite difference method Eq. (fT3]) and Horowitz-Hubeny method 
Eq.()23l). The vertical line separates the stable and unstable regions. 

In the stable region, the behavior of the real part of the quasinormal frequencies as 
L increases is quite similar in both approaches. However, after the phase transition, the 
Horowitz-Hubeny method gives a rather different result. A possible explanation for this 
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Figure 5: Evolution of ^(t) with isolated second mode (left) and second mode behaviours for 
several L (right), (color online) 
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Figure 6: Second mode frequencies (color online). 



behavior could be due to the fact that in the neighborhood of the phase transition point L*, 
the real part of the frequencies computed through this method may be the frequencies of 
the second mode. Indeed, zero frequencies do not show up in the Horowitz-Hubeny method 
due to the its limitations as already discussed in 15] . We also found some frequencies with 
positive imaginary part using the Horowitz-Hubeny method. However, it does not seem to 
converge as we increase the truncation order N. Therefore, we cannot rely on the results in 
this region. 
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Figure 7: Comparison between real frequencies (left) and imaginary frequencies (right), (color 
online) 



According to [19|, polynomials with large terms (e.g. ~ 10 18 ) can be pathologic. A 
vanishingly small change in a very small term can drastically change the roots in such cases. 
In the Horowitz-Hubeny method, the larger the truncation order N, the larger the terms 
of the polynomial are. Thus, for sufficiently large values of N, even small numerical errors 
from polynomial calculation can lead to a completely different root. That explains why the 
roots stop converging after a certain N. In this work, the first derivative term in (fTOj) implies 
in larger terms in the polynomial and the algorithm stops converging for N smaller than 
what is expected for systems with q = 0. The approximate positions of the zeros as well as 
stability are confirmed as well. 

Therefore, we conclude that there are exactly two phases in the model, a stable phase for 
small values of L (normal phase) and a unstable superconducting phase for L > 1.28, or for 
temperatures T < 0.21. We expect such results and phase structure to remain true in other 
holographic models (as e.g. with different topologies). 
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